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Absence of Hybridization Gap in Heavy Electron Systems and Analysis of YbAla in 
terms of Nearly Free Electron Conduction Band 
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In the analysis of the heavy electron systems, theoretical models with c-f hybridization gap 
are often used. We point out that such a gap does not exist and the simple picture with 
the hybridization gap is misleading in the metallic systems, and present a correct picture by 
explicitly constructing an effective band model of Yb AI3 . Hamiltonian consists of a nearly free 
electron model for conduction bands which hybridize with localized f-electrons, and includes 
only a few parameters. Density of states, Sommerfeld coefficient, f-electron number and optical 
conductivity are calculated and compared with the band calculations and the experiments. 
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1. Introduction 

The 4f-electron systems in rare-earths have attracted 
much interest since the duality between the itinerant and 
the locaHzed properties of f-electrons gives rise to rich 
variety of phenomena.^ In particular, Ce and Yb com- 
pounds have been extensively studied for several decades 
because Ce has only one 4f-electron and Yb only one 4f- 
hole. The periodic Anderson model (PAM) is often em- 
ployed in order to analyze the properties of these com- 
pounds. The most simplified PAM neglects the orbital 
degeneracies of both conduction and f-electron states, 
and is given by 



kcr 



kcr 



(1) 



where cl^^ (/ko-) ^ creation operator of a conduction 
(f) electron with the wave vector k and the spin a. 
represents the local f-electron number operator {=fl^ficr) 
J> . with the site index i. is the conduction band energy 



and Ef the f-electron energy level. Vko- represents the hy- 
bridization between the conduction and f-bands (c-f hy- 
bridization) and U is the Coulomb interaction between 
the f-electrons. kcr-dependencc of Vka is often neglected. 
This Hamiltonian is regarded as the minimal model to 
describe the 4f- and 5f-electron systems and has been 
partially successful in qualitatively capturing the physi- 
cal properties. Especially, the heavy mass is explained by 
the combination of the position of the Fermi level inter- 
secting the nearly flat part of the hybridized bands (see 
Fig. 1(a)) and the renormalization due to strong corre- 
lations in f-electrons. Furthermore, a hybridization gap 
appears and the width is of the order of V for the non- 
interacting case. Using this simple model, Okamura, et 
al. found a universal scaling property on the peak posi- 
tion and the hybridization strength over many Ce and 
Yb compounds.^ 

The degeneracy of the f-states is very important. To 
incorporate it, the following model has been investigated 



(a) PAM 




(b) PAM with orbital degeneracy 



Fig. 1. Schematic band dispersions (a) usual PAM and (b) PAM 
with orbital degeneracy of f-electrons. Dashed line stands for the 
Fermi level. 



in detail:'^ 



^ = y"4cLck. + Ei?f/L/, 
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(2) 



kcr A/ 



i,M>M' 



Here, the conduction electrons have been treated as free 
plane waves, and M denotes the z-component of the to- 
tal angular momentum J of the f-electrons. Anderson, 
however, pointed out that only two of (2 J -I- l)-fold de- 
generate f-states have finite mixing matrix elements (see 
eq. (21) later) with the plane waves, so that {2J+ 1) — 2 
f-states remain unhybridized and stay in the gap (see 
Fig. 1(b)). ^ Using this model, Kontani derived the grand 
Kadowaki- Woods relation,^ which improves the original 
Kadowaki- Woods relation by taking into account the ef- 
fect of the degeneracy of the f-orbitals. Comparison with 
the experiments is made in detail in ref. 6. In this the- 
ory, the hybridization was approximated as a constant 
independent of k or M. 

The k-dependence of Vko-M was investigated e.g. by 
Hanzawa et al.J Ikeda and Miyake,^ etc., but the con- 
duction electrons were always treated as free plane waves. 

Completely different approach to heavy electron sys- 
tems is taken by using the tight-binding model to both 
the conduction and f-electrons. For example, the local 
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density approximation (LDA) energy band of the most 
typical Kondo insulator YbBi2 was reproduced very well 
by this approach, and it was clarified that the degenera- 
cies of both the conduction and f-electrons are essentially 
important in the opening of the gap in this material.® 
Anisotropy in the magnetic field-induced insulator-to- 
metal transition was explained^'' and correlation effect 
on the optical conductivity was also investigated.^-'^ Ap- 
plication to Ce-compounds was done by Maehira, et al}^ 
Fermi surfaces of heavy electron systems are calculated 
via LDA and agreement with the de Haas van Alphen 
(dHvA) experiments has been achieved in several com- 
pounds.^'^ LDA calculation, however, does not include 
dynamical many-body effect, which is of much impor- 
tance in many cases. Recently, combination of the LDA 
band calculation and the dynamical mean field theory 
(DMFT) has been rapidly developed.^'* Many rare-earth 
compounds, however, have extremely small energy scale, 
so that it is not easy to carry out the LDA+DMFT cal- 
culations. 

Therefore, it is desirable to construct an effective band 
model for the 4f-electron systems as a starting point to 
include strong correlations. It is already realized by using 
the tight-binding model as mentioned above, but there 
exists other type of heavy electron material which can 
be more suitably described by the nearly free electron 
(NFE) conduction band than the tight-binding model. 
The LDA result indicates that the conduction band 
structure of the heavy electron metal YbAla has the NFE 
character, which is shown in Fig. 2.^^ The band disper- 
sions show the apparent parabolic structure in the low 
energy range below the Fermi energy E-p. 

In this paper, we focus on the heavy electron metal 
YbAla with the NFE conduction band and construct an 
effective Hamiltonian, which may become a good start- 
ing point to systematically discuss correlation effects on 
the physical quantities in the future. To take account of 
the effect of a weak periodic potential for the conduc- 
tion electrons, we use the NFE scheme with the empty- 
core pseudo-potential method,^^'^'' which has been suc- 
cessfully applied to simple metals and semiconductors. 
Furthermore, we accurately consider the hybridization 
effect between the NFE conduction bands and the local- 
ized f-bands, which considerably modifies the previous 
studies. Using this Hamiltonian, we calculate the band 
dispersions, the density of states (DOS), the Sommerfeld 
coefficient and the optical conductivity. We discuss the 
physical properties of typical rare-earth heavy fermion 
metal YbAls by comparing with the LDA and experi- 
mental results. 

YbAla has CusAu-type crystal structure and belongs 
to the simple cubic structure as the Bravais lattice with a 
lattice constant 4.2 A. The effective mass observed from 
the dHvA measurement"'^^ is 14 ~ 24 mg where mg is the 
rest mass of an electron. The Sommerfeld coefficient is 
7 ~ 40 mJ • K~^-mol~"'^."'®'^° The mean valence is in the 
range from 2.65 to 2.8,^^'^^ which indicates that YbAla 
is a valence-fluctuation compound. Note that ref. 22 in- 
dicated the limitation of the application of the single 
impurity Anderson model to the periodic systems in the 
analysis of photoemission spectra. The single site Kondo 
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Fig. 2. LDA band structure of YbAls.^^ Upper panel (a) and 
center panel (b) show the overall band structure and its mag- 
nification near the Fermi level Ep, respectively. Lower panel (c) 
shows the result with the upward-shift of the 4f-levels by 0.2 
Ry.i« 



temperature Tk is 500 ^ 700 K.^® The coherence tem- 
perature is about Tcoh ~ 35 K, below which the electric 
resistivity follows the relation: p cx 19,23 ^j^jg result 
indicates explicitly that YbAla exhibits the Fermi liq- 
uid behavior below Tcoh- The optical conductivity shows 
the pseudo-gap structure at low energy region and the 
mid-infrared peak at ftw ~ 0.25 eV.^'* 

Note that the additional shift of the 4f-levels (0.2 Ry) 
must be introduced in the LDA calculation (Fig. 2(c)) 
in order to reproduce the angular dependence of the 
dHvA frequency.^* In comparison with Fig. 2(b), the 
band around F point near the Fermi level shifts to higher 
energy region, and the electron-like Fermi surface disap- 
pears. The bands around M point near the Fermi level 
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shift down, so that the multiply-connected Fermi surface 
appears. In addition to the energy dispersion, the DOS 
has been obtained by LDA.^^ We discuss the accuracy of 
the constructed model by comparing with these results. 

This paper is organized as follows. We construct the 
effective Hamiltonian using the NFE scheme with the 
empty-core pseudo-potential method for the conduction 
bands, and introduce the hybridization with the 4f-states 
in the next section. In Sec. 3, we show the band disper- 
sion, the DOS and the optical conductivity. Summary 
and discussions are given in Sec. 4. 

2. Construction of Effective model 

In this section, we construct the effective Hamilto- 
nian for the band energy of the rare-earth compound, 
YbAls, which consists of three Al and one Yb atoms per 
unit cell. The electronic configuration of each atom is 
[Ne](3s)2(3p)i for Al and [Xe]{Af)^^{5df{6sf for Yb. 
Except the closed shell parts [Ne] and [Xe], there are 
25 electrons per unit cell of YbAls. While it is natural 
to assume that Al ion becomes trivalent, the valence of 
each Yb ion can be divalent or trivalent, where the mean 
configuration of f-electrons is in the range from /^^ to 
/^^. It means that YbAls is a typical valence fluctuation 
compound and shows a metallic state. 

The conduction electrons in YbAls consist of three va- 
lence electrons from each Al atom and 2 ~ 3 valence elec- 
trons from Yb atom, which spread over in the crystal. 
Although the 4f-orbitals of Yb atom have the strongly 
localized character, the hybridization with the conduc- 
tion electrons makes the f-electrons partially itinerant. 
Therefore, the eS'ective Hamiltonian can be written as 
follows 



/band 



(3) 

where y}'^^'^ represents the band energy. H.'^ (H^) is the 
conduction (localized f) electron term and Ti.^^^^ is the 
hybridization between the conduction electrons and the 
f-electrons. Hint represents the Coulomb interaction be- 
tween the f-clc!ctrons. We focus on the coiistruc;tion of 
the effective Hamiltonian for the band energy 7^'^^'^'^ in 
the present study. 

Let us start with the construction of the effective 
Hamiltonian. First, we discuss the conduction electron 
term H'^. The wavefunction of conduction electrons gen- 
erally spreads over on the large region in comparison with 
a lattice constant. Therefore, the Hamiltonian of the con- 
duction electrons is given by 



H'^-—W^ + Vir), 
zm 



(5) 



where ^(r) stands for the periodic potential term due to 
the ions on the lattice sites, which is written as 



(6) 



where v^^ and v^^ are atomic potentials in the Yb and 
Al sites. T represents the coordinate of each Al ion mea- 
sured from the Yb site in the unit cell, and R is the 



translation vector. We neglect effect of Coulomb interac- 
tions between the conduction electrons because of strong 
screening effects. 

The NFE method is employed for the conduction elec- 
trons in order to take account of the effect of the weak 
periodic potential, which is given by 



y(r)=^Vbe'«--. 



(7) 



G stands for the reciprocal lattice vector. Vg represents 
the matrix element between the plane waves. 



Vg = (k|y(r)|k- G) = ^J dry(r)e-' 

= J_ / drz;Yb(r)e-'G-'- 
"0 J 

+ ^5^e-'^-^y^drt;AHr)e-" 



G r 



iG r 



(8) 



where f2 is the volume of the system and flo is the volume 
per unit cell. 

Next, let us give the potential explicitly. From the 
LDA band calculation, the conduction bands below Ep 
look like NFE. Above Ep, the contribution from d elec- 
trons deforms the bands from NFE. We particularly fo- 
cus on semi-quantitative description of the conduction 
bands around and below Ep. We apply the empty-core 
pseudo-potential method-*^^' to the description of the 
weak atomic potential in each site, which is given by 



t;"(r) = 







(|r| < r^) 
(|r| > r^), 



(9) 



where r" represents the core radius and Z" is the effec- 
tive valence number (a=Yb, Al), which are adjustable 

parameters. 

However, it is well known that aluminium is a triva- 
lent metallic system and its band dispersion can be de- 
scribed well by the pseudo-potential technique. Thus, 
r^' and Z-^^ are taken from Harrison's textbooks^^' "'^'^ 
in the present study. Namely, we set r^' = 0.61 A and 
Z-^^ = +3. Compared with the Al case, the estimation of 
Yb parameters is not so unique. It is because Yb ion can 
be divalent {Z^^ = +2) or trivalent {Z^^ = +3). We in- 
vestigated the effects on the physical properties for both 
values, and obtained almost similar results. We hereafter 
employ Z^^ = +2. Furthermore, we also calculated the 
band dispersions for several values of the core radius r^^ 
in the range from 0.5 to 3.0 A. Although there was no 
qualitative difference for all cases, the overall band struc- 
ture for rj^ = 2.0 A well reproduced the LDA result. 
Thus, we employ this value of r^^ for a trial, although 
it is larger than the value in refs. 16 and 17. 

The Fourier transformations of the pseudo-potentials 
are written as 



AnZ^e^ cos(gr^) 



(10) 



where wc define q = |q|. 

The redistribution of the conduction electrons gives 
rise to screening of the charge at large distance. The 
screening effect in the pseudo-potential is considered by 
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Fig. 3. Bare and screened pseudo-potentials of Yb (upper panel) 
and Al (lower panel) for YbAl3. The effective valence numbers 



are Z^b = 
A and r^^ 



+2 and 



+3, and the core radii are ri" 



2.0 



: 0.61 A, respectively. 



using the dielectric function of the free electron gas. 
Then, we obtain the screened potential 



e(q) 



e(q)' 

4776^ 



I^o(^f) 

ft 



(11) 



no(a;), (12) 



where Do{Ef) is the DOS at the Fermi level. We define 
X = q/2kp where kp is the Fermi wave number. The 
values of the free electron gas are applied to Do{Ef) and 
kp. Ho (a;) is the three-dimensional Lindhardt function, 



no(x) = i 



1 



2x 



In 



1 



1-x 



(13) 



Since wc start with the divalent state for Yb (Z^^ = 
+2), the free electron gas consists of nine electrons from 
three Al^+ ions and two electrons from one Yb^"*" ion. 
The potentials in the Fourier space are shown in Fig. 3. 
Both figures exhibit explicitly that the amplitude of each 
potential is strongly reduced by the screening effect, in 
particular, in the small |q| region. 

Finally, the pseudo-potential for the conduction elec- 
trons in the NFE method for the whole lattice is given 
as follows 



Gt 



(14) 



Thus, the Hamiltonian for the conduction electrons with 



spin a is written as 

/ eki + Vo 



E 

k 



V, 



Gi — G2 



^^Gi-G2 

ek2 + ^0 



V ^Gi-G„ ^G2-Gn 

Gj and kj with an integer set i 
as 



Vbi-G„ \ 
''^G2-G„ 



• (15) 



ek„+Vb / 
(*i)i2)*3) are defined 



Gi = hGl+i2Gl+hGl, (16) 

ki = k + Gi, (17) 

where Gq {h = 1,2,3) denotes the primitive translation 
vector of the reciprocal lattice. Note here that is in- 
dependent of the spin index a. 

Because of the strong spin-orbit coupHng of f- electrons, 
the local f-electron states are specified by the total an- 
gular momentum J and its z-component M. For Yb ion, 
the local f-electron states consist of J = 7/2 and J = 5/2 
states. While the former has the eight-fold degeneracy 
and is located close to the Fermi level, the latter has the 
six-fold degeneracy and is far from the Fermi level for 
YbAla.^^'^^ The energy difference between the J = 1/2 
and J = 5/2 states is defined as A, which is of the order 
of 0.1 Ry.i'''^''' Although the J = 5/2 state is almost ir- 
relevant to the low-energy physics, we include this state 
in Ti} . The localized f term of the Hamiltonian is written 
as 



= E ^ifiJMfiJM, 



(18) 



i.J.M 



where fjjj^f {Jum) is a creation (annihilation) operator 
of a localized f-electron with the site index i, and the z- 
component M of the total angular momentum J = 7/2 
and 5/2. The 4f-energy level Ef is a tuning parameter in 
the present scheme. Note here that although it is easy to 
take account of the crystal field effects it is not discussed 
in the present study, for simplicity. 

Next, let us consider the matrix element of the hy- 
bridization between the conduction electron and the f- 
electron. The plane wave and the localized f-electron can 
be expanded by the spherical harmonics, respectively. 
Considering the orthogonality of the spherical harmon- 
ics, the matrix element of the hybridization term is writ- 
ten as 

VuajM = (ka|H^>"^|/kjM) 

(''Imcr^lm i&k , f^k) Vfc„i , 



Vfen 



47r(-l)' 



(19) 



dr r'^Rniir)v{r)ji{kr),{20) 



where a stands for the spin index of the conduction 
electron. |/k jm) denotes the Fourier transformation of 
I/iJm)- Yim{(^k: (pk) is the spherical harmonics with the 
orbital angular momentum / and its z-component m 
where and (p^ denote the direction of the wave vector 
k. ji{kr) is the spherical Bcssel function.^' For simplic- 
ity, we regard Vkni as a constant parameter V. Although 
Vkcr JM is independent of the radial part of the wave vec- 
tor k, the angular dependence of k still remains, af^^ is 
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the Clcbsch-Gordan coefficient. For the f-electrons with 
I = 3, J = 7/2 and 5/2, it is given by 



(a) 



C-lma — 



-a,J\{l-Ma)5m,M-^ (j = |) , 



(21) 



where a stands for 1 (| spin) and — 1 (i spin), respec- 
tively. 

Therefore, substituting eq. (21) for eq. (19), VkcrjM is 
written as^ 



JM- 



i(| + Ma)y3,M-f(^k,0k)V^ (J=i) 

a^i(|-M<7)y3,M-f(^k,0k)^ {J=l)- 

Finally, the effective Hamiltonian for the band energy 
is written as 



n 



band 









w 



-^hyb 



(23) 



Summarizing this section, the effective Hamiltonian for 
the band energy is constructed by the combination of 
the NFE scheme for the conduction electrons and the 
hybridization with the localized f-electrons. This Hamil- 
tonian includes only a few parameters (the hybridization 
V and the 4f- level £"/ in addition to and r"), which 
can be estimated from the LDA and experimental results, 
as will be is discussed in the next section. 



3. Results 

Although all reciprocal lattice vectors should be in- 
cluded in eq. (15) in the NFE method, the present calcu- 
lations actually converge with i^=— 3 to -|-3 (/x = 1, 2, 3) 
in eqs. (15) and (16). Then, the size of the Hamiltonian 
matrix Ti)^^'^'^ becomes 700 (= 7^ x 2 + 14) for each wave 
vector k. Diagonalizing this matrix numerically, we can 
obtain all eigenvalues for each k. Although we divide 
the first Brillouin zone into N ~ 120'^ pieces to perform 
k-integration, the calculated physical quantities almost 
converge at A'' = 60^. The actual calculations are carried 
out within 1/48 of the first Brillouin Zone. In summary, 
we use the following parameters in the present study, the 
energy difference between the J = 7/2 and J = 5/2 levels 



arc chosen as A = i?. 



7/2 



5/2 



0.1 Ry. The pseudo- 



potential parameters are r^^ = 0.61 A and r^'° = 2.0 A. 
The valences of the ions are simply assumed as Z^^ = +2 
and = +3, respectively. The Yb 4f J = 7/2 level 

7 /2 

Ef and the hybridization amplitude V are the main 
tuning parameters. Except the result of the optical con- 
ductivity, we discuss the physical properties at absolute 
zero temperature. 

Figure 4 shows the band dispersions of YbAla obtained 

7 /2 

by the present scheme for E^' = 0.59 Ry and V = 0.03 
Ry. Note here that the overall band dispersion is shifted 
by 0.48 Ry in order to adjust the origin of our band to 
that of the LDA band (Fig. 2). The almost flat bands 
at ^ 0.59 Ry and at ^ 0.49 Ry correspond to the Yb 
4f J = 7/2 and J — 5/2 levels, respectively. However, 
because of the hybridization term, the 4f-bands become 
dispersive, which is shown in the lower panel of Fig. 4 for 



OS 



c 
W 



0.5 





V 
































3 



M R 



(22) 




M X 



M R 



Fig. 4. 

-,7/2 



Energy band dispersions of YbAls for V = 0.03 Ry and 



Ey^ = 0.59 Ry. The horizontal dashed-line stands for the Fermi 
level. The upper panel (a) shows the overall structure and the 
lower panel (b) is the extreme magnification close to the Fermi 
level (Bp ~ 0.599 Ry). 



the J = 7/2 band. The Fermi level is located at a little 
above the J — 7/2 level and intersects some dispersive 
bands, which explicitly indicates that the obtained result 
describes the metallic state of YbAls. Furthermore, we 
stress that the hybridization gap does not appear within 
the present treatment of the c-f hybridization and the 
parameter sets. 

It seems that although the details are different from 
that obtained by LDA, the overall structure around and 
below the Fermi level corresponds well to the LDA result 
(Fig. 2(a)). The physical properties mainly depend on 
the low-energy electronic structure. It is consistent with 
the LDA result (Fig. 2(b)) except that close to M point. 
There are small closed hole-like Fermi surfaces on F-R 
line and closed electron-like Fermi surfaces are located 
around F and R points, which are consistent with those 
of the LDA band. However, the neck at M point in the 
LDA band (Fig. 2(b)) is not reproduced. 

The band structures are calculated with other values of 
V in order to consider the robustness of the Fermi surface 
shapes. Although the J = 7/2 band becomes dispersive 
with an increase of the amplitude of the hybridization, 
it seems that there is no topological change of the Fermi 

7 /2 

surface for E^ < 0.64 Ry. It is because the main J = 
7/2 band does not directly intersect the Fermi level. 
With an increase of E^"^ level to higher energy region, 

7/2 

it is approaching to the Fermi level. For E^ > 0.64 Ry, 
the electron-like Fermi surface around F point vanishes, 
which is shown in Fig. 5. This band structure is in partial 
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Fig. 5. Band structure around the Fermi level (dashed line) for 
V = 0.03 Ry and E^'^ = 0.64 Ry. 



agreement with the LDA band with the upward-shift of 
the 4f- levels (Fig. 2(c)). Note that the hybridization gap 
does not appear in this case, either. 
The number of f-electrons is given by 

. 1 ^ 



N ^ ' 

kaJM 



,JM\1 



(24) 



where E^a and wjj^ are the eigenvalue and eigenvec- 
tor of Ti,^^'^'^ with the wave vector k and band index a. 
The results for several choices of the hybridization V are 
shown in Fig. 6. With an increase of the J = 7/2 level 
eJ^"^ and/or the hybridization V, the deviation from f^* 
configuration becomes conspicuous. Compared with the 
experimental mean- values of Yb valence number (2.65 

~ 2.8), eJ^^ is expected as 0.59 ~ 0.60 Ry. While the 
band structure in Fig. 5 can partially reproduce the LDA 
result with the upward-shift of the 4f-levels (Fig. 2(c)), 
the number of f-electrons becomes too small (nf ~ 12.72 
< 13). Thus, we hereafter focus on the parameter set in 
Fig. 4 (rif - 13.57 for eJ^^ = 0.59 Ry and V = 0.03 Ry). 

Next, let us study the property of the DOS and the 
number of states (NOS). Figure 7 shows the total DOS, 
the partial f-DOS and the NOS, which are defined as 

^*°*(^) = ^E^('^-^ka), (25) 



13.5 



■i 



<^ V=0.02 
*- V=0.G3 
^ V=0.G4 
V=0.05 



12.5 




Ef''' (Ry) 



Fig. 6. Number of f-electrons as a function of J ■ 
for several choices of hybridization V. 
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Fig. 7. (a) DOS for V = 0.03 Ry and eJ^^ = 0.59 Ry. The solid, 
dashed and dotted lines stand for total, J = 7/2 and ,/ = 5/2 
DOS's, respectively. The inset shows the magnification around 
the J = 7/2 level. The panel (b) shows NOS where the dashed 
line stands for the total particle number of YbAls per unit cell, 
which corresponds to 25. The dash-dotted line represents the 
Fermi level in both panels. 



^J'M = ^EK'l''^(^-^ka), (26) 
kaM 

N{uj) = r dw'£)*°*(a;')- (27) 

J —CO 

The overall structure of DOS is roughly similar to the 
result of ref. 25. Two sharp peaks appear around the 
J = 7/2 and J = 5/2 levels because of the strongly lo- 
calized character of the f-electrons. However, both peaks 
are broadened due to the c-f hybridization. With an in- 
crease of the hybridization amplitude, the peak around 
the J = 7/2 level splits into two parts, which is shown 
in the inset of Fig. 7(a). The energy difference between 
those peaks is about 0.005 Ry for the given parameters. 
It is important that the hybridization gap does not ex- 
plicitly appear because of the angular dependence of k 
in the c-f hybridization, the periodic potential and the 
band-folding into the first Brillouin zone. Of course, al- 
though the further increase of the hybridization ampli- 
tude gives rise to the growth of the pseudo-gap, an ex- 
plicit gap never appears within the range of the investi- 
gated parameters. Since the partial 4f-D0S has the tail 
due to the c-f hybridization, the Fermi level partially 
intersects the J = 7/2 bands. Thus, the number of f- 
electrons per unit cell deviates from integer (n^ = 14) 
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Fig. 8. Optical conductivity for a few choices of J = 7/2 level and 
V = 0.03 Ry at absolute zero temperature. The dash-dotted line 

7/2 

stands for the result for E^' — + — oo and V = 0. Drude parts 
are not shown. The inset shows the energy difference between 
the Fermi level E-p and the J = 7/2 level E^^ . 

except the V = Q case. Figure 7(b) shows the NOS for 
reference, in which the jumps appear around the J = 7/2 
and J = 5/2 levels. The presence of the c-f hybridization 
smears the jumps. Similar to the band dispersion, the 
result of the DOS also indicates that YbAls is in a va- 
lence fluctuating state, which is in agreement with the 
experimental results. 

The Sommerfeld coefRcient 7 is written as 

7=yfci^*°'(£;F)iVA, (28) 

where A^a is Avogadro's number. For E^' = 0.59 Ry 
and V = 0.03 Ry, 7 is about 8.55 mJ- R-^-mor^ For 
other values of the parameters, 7 is in the range from 
5 to 10 mJ- K~^-mol~^ . These values are smaller than 
the experimental value (~ 40 mJ- K~^-mol~^) since the 
effect of electron correlations between the f-electrons is 
not taken into account, which will be discussed later. 

Finally, let us investigate the optical conductivity (j{uj) 
and its temperature dependence. By using Kubo formula, 
the contribution of the direct transition, except Drude 
part, is given by 

Rea"(a;) = -^ ^ {3f{\.)f{f{E^^)-f{E^p)} 

xS{u;-{E^f3-Ek^)), (29) 

where j"'^(k) stands for the current matrix element be- 
tween the band index a and 0. f{E) is the Fermi dis- 
tribution function. For simplicity, we neglect the wave 
vector and band index dependences of the current term 

in the present study. Namely, we set j"^{]s.) = constant. 

Figure 8 shows the conductivity for a few choices of 
7/2 

Ef at absolute zero temperature, in which Drude parts 

are not included. To consider contributions of the f- 

7/2 

electrons, we also calculate the conductivity for E^ 

7/2 

—00 and V = 0. Except the Ef ^ — 00 case, the peaks 
appear at low energy region (w ~ 0.01 to 0.03 Ry, which 

corresponds to the range from 0.15 to 0.4 cV). With a de- 

7 /2 

crease of the E^ level, the peak structures are smeared 



and the peak positions shift to the high energy region. 

7/2 

For Ef —00, the amplitude is extremely smaller 

than those of other values of eJ^^. Thus, we can con- 
clude that these peaks at w ~ 0.01 to 0.03 Ry are mainly 
due to the excitations from the J = 7/2 bands to the con- 
duction bands above E-p with the same k, in which the 
peak position for E^"^ = 0.59 Ry corresponds roughly to 
the energy difference between the center of gravity posi- 
tion of the J = 7/2 DOS and the Fermi level (shown in 
the inset of Fig. 7(a)). The peak due to the excitations 
from the J = 5/2 bands appears at w ~ 0.1 Ry. Note 
that we replace the delta function with a gaussian func- 
tion in eq. (29) on the numerical calculations, so that 
this broadening effect gives rise to the sharp peaks close 
to w = 0. With a decrease of a half width of the gaussian 
function, the peak amplitude becomes small. 

The optical c;onductivity observed in the experiment 
shows the pseudo-gap structure at very low energy region 
and the mid-infrared peak around tuv = 0.25 eV.^^ The 
origin of the mid-infrared peak has been attributed to the 
excitations across the c-f hybridization gap, seen in the 
simple PAM.^^ However, in our model, the hybridization 
gap does not exist (shown in Fig. 4 and the inset of Fig. 
7(a)). 

The structure of the calculated conductivity is consis- 
tent with that of experimental result. The peaks found 
in the range from 0.01 to 0.03 Ry seem to correspond to 
the mid-infrared ones observed in the experiment. The 
energy difference between the split peaks in the J — 7/2 
DOS (the inset of Fig. 7) is at most 0.005 Ry (~ 0.07 
eV), which is too small for the energy scale of the ob- 
served mid-infrared peak. Of course, with an increase of 
the hybridization amplitude, the energy difference be- 
tween the split peaks in the J = 7/2 DOS is enlarged 
and the tail in the partial f-DOS grows. In this case, 
since the number of f-electrons has to be in the range 
from 13 to 14, the Fermi level shifts to the higher energy 
region away from the J = 7/2 level (shown in the inset 
of Fig. 8). Therefore, the peak position of the conduc- 
tivity also shifts to the higher energy region. Thus, we 
conclude that the origin of the mid-infrared peak is not 
the excitations between the split bands due to the c-f hy- 
bridization but the excitations from the J = 7/2 bands 
to the conduction bands above Ep with the same k. 

The temperature dependence of the conductivity is 
shown in Fig. 9. Drude parts are not included. The 
pseudo-gap structure at low energy region is gradually 
filled with an increase of temperature. The characteris- 
tic temperature for disappearance of the pseudo-gap is 
about 500 K for the present parameters. If correlation ef- 
fects between the f-electrons are included, a Drude peak 
with a finite width will be added at finite temperatures, 
and the width will become larger with an increase of 
temperature. Filling of the pseudo-gap will become faster 
than that of the non-interacting case since the DOS also 
depends on temperature.^^ 

4. Summary and Discussions 

We have constructed the effective model of the band 
energy for heavy fermion compound YbAls. In order to 
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Fig. 9. Temperature dependence of optical conductivity for V = 
0.03 Ry and eJ^^ = 0.59 Ry. 

take account of the effect of the periodic potential reflect- 
ing the crystal lattice structure, the NFE method with 
the empty-core pseudo-potential is employed for the con- 
duction electrons, in which the screening effect is consid- 
ered by using the dielectric function of the free electron 
gas. The localized bases are employed for the 4f-electrons 
because of the strongly localized character. Furthermore, 
the hybridization between the NFE conduction band and 
the f-band is accurately taken into account. Adjusting 
only a few parameters and diagonalizing the obtained 
Hamiltonian, we have calculated the band dispersions, 
the DOS, the Sommerfeld coefficient and the optical con- 
ductivity. 

While the often-used simple PAM and its extensions 
usually yield a hybridization gap, our realistic model ex- 
plicitly denies the existence of the hybridization gap in the 
metallic system. The number of f-electrons deviates from 
integer and covers the experimentally observed range. 
Thus, the valence fluctuating state is reproduced. The 
obtained band structure and the DOS are mostly con- 
sistent with the LDA results. Furthermore, we find that 
the origin of the experimentally observed mid-infrared 
peak is due to the excitations from the dispersive 4f 
J = 7/2 bands to the conduction bands above with 
the same wave vector k, and not the transition across the 
hybridization gap. We conclude that the present effec- 
tive model with only a few parameters may describe the 
physical properties of YbAls. However, by varying the 
parameters, we can only partially reproduce the LDA 
band with the upward-shift of the 4f-levels which is in 
agreement with the dHvA result. In addition, the num- 
ber of 4f-electrons becomes unfortunately too small. This 
discrepancy in the present approach is to be improved in 
a future. Nevertheless, the present approach seems to be 
applicable to a wide range of heavy electron metals, and 
our conclusion that the hybridization gap is absent would 
be true for most cases. 

Finally, we have to discuss correlation effects of the 4f- 
electrons. Although it is well known that the rare-earth 
compounds arc strongly correlated electron systems, we 
focus on the construction of the effective Hamiltonian 
of the band energy in the present study. Correlation ef- 
fects should be taken into account in order to estimate 



physical quantities accurately, which is now under inves- 
tigation. However, the Coulomb interactions are strongly 
screened at low energy region in heavy fermion systems, 
so that the Fermi liquid behavior is often realized at low 
temperature region. For T < 35 K, YbAls shows the 
Fermi liquid property. Thus, if the 4f- levels and the c-f 
hybridization amplitude employed in the present study 
are regarded as the renormalized quantities, correlation 
effects may be taken into account partially. Of course, 
the physical properties caused by dynamical quantum 
fluctuations are not included. Our model may be a good 
starting point to systematically consider such correlation 
effects. 
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